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A Sinple  Batch  Epidemic  Process 
by 

L.  Billard,  H.  Lacayo  B N.  A.  Langberg 
Florida  State  University 

SUKMARY 

A simplifying  assuiq)tion  for  an  epidemic  process  is  that  at  most  one  person 
may  becooM  infected  at  any  one  time.  However,  it  is  quite  conceivable  that  when  an 
infected  person  makes  simultaneous  contact  with  two  people,  both  could  become  in- 
fected. In  this  note  we  introduce  the  concept  of  a batch  epidemic  process  in  which 
infection  can  occur  in  batches  of  one  or  two.  The  distribution  of  the  number  of 
infectives  present  at  any  time  is  derived.  The  concept  can  be  extended  to  in- 
clude batches  of  any  size. 
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1.  INTRODUCTION 

The  slaqple  epidemic  model  describes  the  process  in  which  a population  of 
fixed  site  N consists  of  susceptibles  and  infectives.  A basic  assum>tion  is 
that  at  any  particular  instant  in  time  at  most  one  susceptible  individual  can 
becoaM  infected.  However,  it  is  quite  conceivable  that  whm  an  infected  person 
makes  simultaneous  contact  with  two  people  both  could  become  infected.  More 
generally,  susceptibles  could  become  infectives  in  batches  of  any  size.  For 
simplicity,  we  shall  consider  here  the  case  in  which  Infection  occurs  in  batches 
of  one  or  two  with  probability  a and  (1  - a)  respectively. 

Let  X(t)  denote  the  number  of  infectives  at  tine  t.  The  infinitesimal 
transition  rates  for  our  model  are  given  by 

Pr{X(t  ♦ h)  ■ X ♦ l)X(t)  ■ x)  ■ ow^h  ♦ o(h) 


Pr{X(t  ♦ h)  ■ X ♦ 2|X(t)  ■ x)  ■ (1  - Q)Mjjb  ♦ o(b) 


and 


Pr{X(t  ♦ h)  ■ x|X(t)  ■ x)  ■ 1 - * o(h) 


(1) 


where  lij^,  x ■ 1,  ...,  N,  are  arbitrary  but  distinct  infection  rates.  In  the 
classical  simple  epidemic  model,  a s 1 and  0x(N  - x)  where  6 is  the  rate 
of  infection.  In  this  work  we  find  expressions  for  the  probability  distribution 
of  X(t)  for  the  above  simple  batch  epidemic  process.  Note  that  the  underlying 
model,  from  (1),  satisfies  the  Chapckan-’Co>icgorov  equations 

- -M„p„(t)  ♦ aWn.iP^^.lCt)  ♦ (1  - «)»«n.2Pn-2^*^*  (2) 

for  n ■ I N with  w ■ 0 for  n 2 N or  n < 1 ard  wliere 

n 

/ 

p^(t)  - Pr{X(t)  - n) 

*^n  t 
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2.  SOME  PRELIMINARIES 

I 

Recently,  Billar<i^  Lacayo  and  Langberg  (1979)  found  a compact  derivation 
for  the  classical  simple  epidemic  process  by  focusing  on  the  interinfection 
time  between  successive  occurrences  of  infectives.  The  probability  distribu- 
tion for  the  number  of  infectives  present  at  any  time  t was  expressed  in  terms 
of  a linear  combination  of  exponential  terms  with  parameters  where  is 

defined  as  in  Section  1.  Since  in  the  classical  model  infectives  occur  in 

batches  of  one  only,  the  patn  followed  by  successively  occuring  values  of  x is 
simply  1,  2,  ....  N.  However,  when  infectives  can  occur  in  batches  of  one  or 

two,  there  are  many  possible  paths  that  can  be  followed.  Accordingly,  we  intro- 
duce some  notation  to  describe  the  path  possibilities. 

Let  Uj,  ...,  he  distinct  positive  real  numbers.  Define  M(n,  i,  j)  as 
the  set  of  all  subsets  of  with  j of  them  deleted  but  such 

that  no  two  deleted  elements  have  consecutive  subscripts.  Thus,  for  example, 

M(6,  1,  2)  - {(U-,  w.),  (w,,  w.),  (w,,  M-)}.  Let  A » (u  , , .. .,  Uy  ) 

^ ^ ^ ^ ^ ^1  V2  %-i-j-l 

denote  a typical  element  of  M(n,  i,  j).  In  essence,  if  we  are  given  i infec- 
tives, the  subscripts  of  the  components  of  A represent  a particular  path  fol- 
lowed to  reach  the  state  corresponding  to  n infectives.  Since  the  end  points, 
i and  n,  are  always  present,  these  subscripts  are  suppressed.  For  example. 
Figure  1 indicates  the  possible  paths  from  i > 1 to  n > 8 whenever  j ■ 3 jumps 


of  size  2 occur. 
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i ■ 1 n ■ 8 

»'4*  “bJ 
•'4*  ‘'6^ 

“5* 

“5* 

FIGURE  1 

Paths  for  M(8,  1,  3) 

Note  that,  in  general,  j - 0,  1,  ...,  [(n>i)/23  where  [x]  is  the  greatest  inte> 
ger  less  than  or  equal  to  x.  Let  us  further  define  M(n,  i,  j;  to  be 

those  eleaents  A of  M(n,  1,  j)  such  that  and  let  M(n,  i,  j; 

siadlarly  defined.  Thus,  in  our  exaaqile,  M(8,  1,  3;  equals  by)} 


A - (Mj* 


T^>v 


A - (Pj. 


A ■ (W.. 


A - (Wj. 


since  only  one  pt'.th  of  M(8,  1,3)  contains  Py.  Let 


♦ (n.  i,  j)  - I 

A*M(n,i,J)  V'' 


Pi) 


(3) 


Finally,  let  4(n,  i,  j;  be  the  sun  of  those  terms  of  4(n,  1,  j)  that 

have  the  factor  of  maxlnal  subscript,  with  4(n,  i,  j;  b,j.2)  sinilarly 

defined. 


For  distinct  positive  real  nun^ers  p^,  ...,  p^, 

♦(n  ♦ 1,  i,  J)  ■ p^(pj^  - pj)"'^(n,  i,  j)  ♦ w„_i(Pn  - P^)”^4(n  - 1,  i,  j - 1)  (4) 

where  ^(n,  i,  j)  is  defined  as  in  (3). 

Proof. 

Since  by  definition  the  sets  A of  M(n  ♦ 1,  i,  j)  cannot  have  both  p^  and 


p^_j  deleted,  it  follows  that 


I 
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5.  A SIMPLE  BATQi  EPIDEMIC 


L0t  us  now  presont  the  probability  distribution  of  X(t)  the  nuaber  of  infectives. 


Theorea  1 


Consider  a siaple  batch  epideaie  in  a closed  population  of  sise  N with  transi- 
tion probabilities  as  given  in  (1).  Assuae  the  ^laraiMtors  ■ 1*  •••>  N are 

distinct  and  that  X(0)  • 1.  Then,  the  probability  distribution  of  the  nuaber  of 
infectives  present  at  tine  t is 


Pn“>  ■ ' ” ■ * " • 


where 


'l.!  * *• 


*n  i • I,  J),  * ; 5 ; - *■ 

n,i  »*n  “ **1  j.o  " **  •••• 


1 


n * 2,  •••,  N, 


and  where  the  fCn,  i,  j)  are  defined  in  (5),  and  where  d^  0 whenever  a > n. 

n ,n 

Proof. 

The  result  for  n ■ 1 is  clearly  true. 

he  prove  the  result  for  n ■ 2,  ...,  N by  induction.  We  first  take  n ■ 2. 
Substituting  n ■ 2 into  equation  (2) , we  obtain 

-Ml*  -u^t 

PjCt)  ■ Caii|/(W2  " Wi)3e  - - Pj)3e 

Equation  (8)  implies  the  result  is  true  for  n ■ 2. 

Now  assuae  the  result  is  true  for  n ■ 3,  ...,  n • 1.  We  wish  to  check  the  case 

n ■ n.  Solving  (2)  for  n > n,  by  using  the  integrating  factor  expCw^t) , and  (7) 

for  r.  - 1 , n - 2 we  obtain 


P„(t)  - i taVlVl.i*  - ‘‘i^* 

i«l  » 


* iJl  ^ ■ “^V2V2.i^^^^n  ■ 


■*'n" 

e , . , 


Thus,  to  prove  the  result,  we  need  only  show  that  the  coefficient  of  exp(-ij^t)  in 
(9)  is  indeed  d . . To  do  this  we  substitute  for  d . . and  d . . using  (8) 

HgX  n^lfX  II**  gX 

which  were  assumed  to  be  true.  Let  us  assume  n - i ■ 2ffl  -t-  1 is  odd.  (The 
proof  is  analogous  for  n - i even) . Thus , on  a proper  rearrangement  of  the  terms , 
the  coefficient  of  exp(-Uj^t)  becomes 


Wid. . r ^ . u , 

i IX  n*x  i\* i A I « 

a ♦(!!  - 1,  Xg  o) 

J*U<  U *11* 

n '^i  ’^n-l 


C(n-i)/2] 


j.n-i-2jj  *^n-l 


♦ (n  - 1,  i,  j) 


Vl  ■ ‘'i  ^-2- 


♦ (n  -2,  i,  j - 


I (1  - a)Ja’ 


But,  from  (6)  and  (3),  we  have 

Wn./C»'n-1  ■ “ ♦f*'*  °*  V 

(10) 

• ♦(n,  i,  o). 

Substituting  (10)  and  using  Lenma  1,  yields  the  result  for  n ■ n,  and  the  theorem  is 
proved. 

Note  that  the  coefficients  d , satisfy  the  following  recursive  relationship 

n ,1 

‘*n,i  “ ^®Vl‘*n-l,i  * - “>V2V2.i^/t»'n  ' »'i^  * " " ^ 

•vith  = o -.aei.ever  m > n.  iTiis  foliov;s  iu.aoaiatcly  vrom  (9). 

With  the  expressions  for  the  probability  distribution,  we  may  now  derive  other 


quantities  of  interest,  such  as  the  mean  number  of  infect! ves,  duration  time  of  the 
epidemic,  the  size  of  the  epidemic,  etc. 
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